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Abstract 

We show that the planar spiral phase of the 2D Hubbard model at low 
doping, x, is unstable towards a noncoplanar spin configuration. The novel 
equilibrium state we found at low doping is incommensurate with the inverse 
pitch of the spiral varying as y/x, but nevertheless has a dominant peak in 
the susceptibility at (tt,7t). Relevance to the NMR and neutron scattering 
experiments in La2- x Sr x CuO± is discussed. 
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Magnetic properties of the Cu02 layers in the high temperature superconductors have 
been recently attracting an intense interest [|1] as magnetism is possibly a major contributor 
to the mechanism of superconductivity. Among cuprate superconductors, Lg^CuC^ holds a 
special place as the most simple and extensively experimentally studied compound. There 
are numerous reasons to believe that magnetic properties of La<i- x Sr x CuO{ are quantita- 
tively captured by the 2D Hubbard model with chiefly the nearest-neighbor hopping 

H=~tJ2 4,a a j,a + U n i,T n i,l (!) 
i,j i 

Here a is a spin index, and n = a^a. At half-filling, the ground state of the 2D Hubbard 
model exhibits a long range commensurate Neel order. However, the introduction of even a 
very small number of holes into the system brings about a dramatic change in the nature of 
the ground state as the T = static Pauli susceptibility associated with vacancies does not 
depend on the carrier concentration in two spatial dimensions. 

Shraiman and Siggia pointed out |2| that holes introduced into a commensurate antifer- 
romagnet give rise to a long-range dipolar distortion of the staggered magnetization. In the 
simplest scenario, this leads to a spiral spin configuration with the momentum (ir,Q). The 
incommensurate (ir, Q) spin-density-wave phase was also obtained in the early perturbative 
studies of the Hubbard model with small U || and in several other mean-field |||| and self- 
consistent H calculations. Experimental situation in La 2 - x Sr : x CuO '4 is however far from 
being clear. On one hand, the inelasic neutron scattering experiments have shown that 
x" is peaked at incommensurate (tt, Q) and symmetry related wave vectors, with Q varying 
as tt — Q pa 2ttx. On the other hand, the longitudinal spin-lattice relaxation rate data for 
17 O |8j , obtained at far smaller frequencies cannot be fitted by the form of x" inferred from 
the neutron scattering experiments. 

In this communication we use the RPA formalism of Schrieffer, Wen and Zhang 
and study in detail the structure of magnetic correlations in the Hubbard model at small 
doping. We will show on the basis of an analysis of low frequency bosonic modes that 
the (tt, Q) spiral state is actually unstable with respect to a more complex noncoplanar 
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spin configuration which, though incommensurate, has a dominant peak in the dynamical 
susceptibility precisely at (ir, ir). 

The analysis presented here is related to other works on incommensurate magnetic phases 
at finite doping. Previous mean-field studies of the Hubbard and t — J models 0,f§|| have 
shown that in a certain range of parameters, the (n, Q) spiral state has lower energy than 
both the Neel state and the spiral state with the pitch in both spatial directions. Our energy 
analysis is consistent with their results. Shraiman and Siggia |I| developed a macroscopic 
theory of the bosonic excitations in the spiral phase. To the lowest order in density, they 
found a peculiar degeneracy in the ground state energy for the planar spiral state and for 
a whole set of noncoplanar magnetic configurations with the plane of the spiral varying 
in space. This degeneracy is also present in our microscopic calculations. However, the 
further assumption of Shraiman and Siggia that the next-order terms in doping concentration 
stabilize the planar spiral state is inconsistent with our microscopic results for the Hubbard 
model. 

As an input for our analysis, we will need an expression for the energy spectrum, E^, of 



a single hole in a quantum antiferromagnet. The mean- field theory gives E^ = — y e k + A 2 , 
where = —2t(cosk x + cosk y ) and A is a gap which separates valence and conduction 
bands [f|,[lT],|9],[l2| . This energy is obviously degenerate along the whole edge of the magnetic 



Brillouin zone \k x ± k y \ — 7r, where Ek has a maximum. However, both perturbative [13 
and variational |14| studies have shown that this degeneracy does not survive the effects of 
quantum fluctuations, and the actual dispersion has a maximum at four points (±7r/2, ±vt/2) 
in the center of each of the edges of the magnetic Brillouin zone. In the neighborhood of 
these points Ek can be presented as Ek = —A + p 2 /2m\\ + p 2 ± /2m±, where near (tt/2,tt/2), 
P\\ = (k x — k y )/2, p± = {k x + k y )/2. Self-consistent calculations predict that at large U both 
masses scale as inverse bandwidth J = At 2 /U, but numerically, m\\ (which is infinite in the 
mean- field theory) is several times larger than m±. 

We now turn to the description of our calculations. Consider first a Neel ordered state 
at small but finite doping. Near Q = (tt, tt), the static transverse susceptibility should 



have a hydro dynamic form |L5] XatiQ ~ Qo) = No/(p s (q — Qo) 2 ), where N is the on- 



site magnetization, and p s is the spin-stiffness. This form of transverse susceptibility is 
reproduced in the RPA formalism by summing up the ladder series of bubble diagrams. At 
half-filling, only bubbles containing valence and conduction fermions are allowed, while at 
finite doping chemical potential moves inside the valence band, and one also has contributions 
to x horn bubbles with only valence fermions. These last contributions are proportional to 
2D Pauli susceptibility, which, as we said, does not depend on carrier concentration. As a 
result, one obtains a finite correction to the spin stiffness already at a very small doping |TT 



Ps = P° s (l-z), z = AT X ^= 2 I^^ (2) 

7T 

Here p° s is the spin stiffness at half filling, and T is the scattering amplitude for two holes. 
The latter is equal to U in the mean-field description, but strong self-energy and vertex 
corrections reduce T to the order of the bandwidth J P JTT|JT6[| , which in turn implies that 



z is simply a dimensionless number. For the rest of the paper, we will consider z as a 
phenomenological input parameter WJ\ . It follows from (E3) that Neel state remains stable at 



finite doping if z < 1, but becomes unstable if z > 1, which implies that we have to consider 
incommensurate spin configurations as possible candidates for the ground states. 



Let us first focus on the two simplest candidates P,[q,[18| : the spiral states with the 
ordering vectors (tv,Q) and (Q,Q). For definiteness, we choose the ordering to be in the 
XY plane such that S§ = S'qCos(QR), and S]^ = 5gsm(QR), where Q is the ordering 
momentum. 

The mean-field analysis for the incommensurate states proceed in the same way as for 
the Neel state. We skip the details of the calculations and focus only on the results. For 
(-7T, Q) state, we obtained Q = n — (U/t)x and 
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E (*,Q) _ = g (i _ z \ (3) 

4 v /mjmj| 

Clearly, the (tt,Q) phase has lower energy than the (tt, 7r) phase for z > 1, exactly where 
the susceptibility of the (ir, tt) state becomes negative. 



For (Q, Q) phase, the inverse pitch Q is the same, but the energy difference is 
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E (Q,Q) _ = (3 _ 2z ) . (4) 

4 v /mlmj[ 

Compairing ([|) and (f|), we observe that the (n,Q) spiral phase has the lowest energy at 
1 < z < 2. This is roughly consistent with the results of other mean-field approaches [@||. 
At z > 2, the (Q, Q) state has the lowest energy. However, we found that 8E^ Q ^ jdx 2 is 
negative at z > 2. This suggests phase separation which is probably unphysical because the 
model neglects long-range Coulomb interaction |2|]. In view of this, we will only consider the 
case 1 < z < 2. 

We now turn to the central topic of our paper, which is the study of collective bosonic ex- 
citations (poles of dynamical susceptibility) in the (ir, Q) spiral phase. The £0(3) symmetry 
of spin rotations is completely broken in the spiral phase, and one, therefore, should expect 
to have three bosonic zero modes with two different spin-wave velocities. The calculations 
of the susceptibilities is straightforward but lengthly because we have to solve a set of four 
coupled Dyson equations. For in-plane (XY) spin fluctuations, a simple symmetry analysis 
predicts that the zero modes are located at q = =)=Q. for and Xq + correspondingly. We 
indeed found that at q = (q+ Q) -C 1, static Xq~ has the form: Xq~ — 2/ (</q 2 (2 — z)), i.e., 
it is positive for z < 2. For the dynamical susceptibility, we found a pole at u = cq with the 



same c (up to 0(x) terms) as at half filling. This is in agreement with other results ||10| , |T8 

We now turn to the magnetic susceptibility Xq Z associated with the fluctuations of the 
plane of spin ordering. This channel is coupled to the charge and in-plane spin fluctuations 
only dynamically, so for the full static susceptibility one has the simple RPA formula x^i 00 = 
0) = x z q z {oJ = 0)/(l — Ux~q Z (u = 0)) where x zz is the bare bubble. From general symmetry 
considerations, we would expect the Goldstone modes in Xq Z to be at q = ±Q. Performing 
calculations to the lowest nontrivial order in the density (i.e., to 0(x 2 )), we have, indeed, 
found that (Xg Z ) _1 at q = ±Q is equal to zero. However, with the same accuracy, we 
also found that the stiffness for excitations near these momentum is equal to zero, which 
in turn means that to order 0(x 2 ), all fluctuation modes between Q and — Q, including 



the mode at (tt, tt), remain degenerate and gapless. A similar degeneracy was found in the 
macroscopic consideration by Shraiman and Siggia ]nj. They further asssumed that 0(x 3 ) 
terms make the spin stiffness at ±Q positive. This scenario yields anomalously small spin- 
wave velocity for out-of-plane fluctuations, v oc sfx, susceptibility at (tt, n) oc a; -3 , and the 
doping-dependent quantum corrections to the order parameter which scale as yfx instead of 
x. 

We, however, calculated the static susceptibility x zz at q = (7r, 7r — q) explicitly to the 
order x 3 and found after tedious algebra: 

where q = tt — Q. We see that while the Goldstone mode at q — ±q survives to order 
0(x 3 ) (as it should), the static susceptibility for out-of-plane fluctuations at and near (tt, tt) 
is negative for z < 2, where mean-field solution favors (tt, Q) phase. This implies that the 
spiral (tt,Q) state is unstable at low doping. Note that all 0(x 3 ) terms come from the 
integration within the hole pockets (the bubbles with conduction and valence fermions yield 
regular corrections in powers of x 2 ). Near the minima, the hole spectrum has a quadratic 
dispersion for an arbitrary form of the hopping term, and we therefore expect eqn. (Q) to 
be valid also for the models with further-neighbor hopping, etc. 

We now address the issue of the true ground state at finite doping. The instability in 
X zz at (vr,7r), implies that the system prefers to have a spontaneous commensurate antifer- 
romagnetic order along Z direction in addition to the incommensurate spin ordering in the 
XY plane. This gives a set of magnetic states which are all noncoplanar (and therefore have 
nonzero chirality), and differ in the ratio between the order parameter amplitudes along Z 
direction and in the XY plane, Fig.|l|. Let us first calculate the ground state energies for this 
set of states. Clearly, we have to introduce two gaps = U (S±), and Aii = U (S\\), where 
(S±) and (S\\) are the magnitudes of the off-plane and in-plane components of the order 
parameter, respectively. Performing the mean-field decoupling and the diagonalization of 
the Hubbard Hamiltonian, we obtained after some algebra 

6 



E A± ^o -Ux 3 



A,, 







)• 









(6) 



where -Ea ± ^o is the ground state energy in the limit when the Z-component of order parame- 
ter tends to zero, and the total order parameter A 2 = Aj_+A| is fixed by the self-consistency 
conditions, A « U/2 [[0|]- It is apparent from (||) that the energy decreases as the ratio 
Aj_/A|| increases. Notice that (i) the energy dependence on the ratio of A is in the order 
0(x 3 ), while to order x 2 , all states from the set are degenerate in energy, and (ii) the r.h.s. 
of (||) contains the same positive factor (2/z) — 1 as the expression for \ zz & t i^,^)- The 
inverse pitch of the spiral is related to x by the requirement that the two self-consistency 
conditions for two order parameters be compatible with each other. We obtained (q = ir — Q) 



U A 

TaT 



1 + 2x 1 



A 2 + AJ 2 
2A 2 



(7) 



As written, Eqn. is valid when q -C Ay/A. Since gAy ~ x (see ([7])), we can keep lowering 
the energy by decreasing Ay as long as Ay ^> \px. When both Ay and q become of the order 
of y/x, AE = E& ± — -Ea x -*o has a more complicated form: 



AE 



Here we introduced A 



Ux 2 



1 + (a 2 + 2)(3 2 - 



f3z 



a 2 + 



3/2 



a 



6 \ \ z, 

aAx 1/2 , q = (3(U/t)x l/2 . Eqn (0) relates a and (3 as: (3 



z((ya 2 + A/z — a)/2. Minimazing now the energy with respect to a we obtain to the lowest 
order in 2 — z: a « Jlz/(12{2 — z)), [3 « 1/a, and 

Ux 2 



AE 



12 / 2^ 2 



(8) 



As expected, AE is negative. Observe also that it scales as x 2 , instead of x 3 as in (||). Clearly 
then, the equilibrium state with Ay ~ ^fx does not belong, strictly speaking, to a set of 
initially degenerate spin configurations, and could be selected already in the calculations 
to order 0(x 2 ). The discovery of the degenerate set of states and of the instability of the 
planar spiral however gave us a hint where to look for the global minimum of the energy. 

Consider next the magnetic susceptibility of the equilibrium state. By virtue of having 
minimum energy it has a positive static susceptibility diverging only at the Goldstone points. 



We also found that noncoplanar ordering gives rise to a strong mixture between in-plane 
and out-of-plane fluctuations, such that the Goldstone modes in xT at q = ±Q give rise 
to a static zero mode in Xq~ & t l 71 "; 71 ")- Another zero mode in x + ~{o) corresponds to the 
spin rotation about the Z axis, and appears at q = — Q, as in the spiral case. Overall, the 
in-plane static susceptibility x + ~ nas two poles 

(q-vr) 2 (q+Q) 2 



where the residues x-n- and xq are proportional to Aj_ and Ajj respectively. Because Ay ~ y/x, 
the residue of the pole at the incommensurate wave vector q = — Q is proportional to the 
hole concentration, and is suppressed with respect to the pole at the commensurate wave 
vector (vr,7r). We also considered dynamical susceptibility and found that the spin wave 
velocity in the zz channel does not contain any smallness in x. However, when z ~ 2, it 
behaves as c ~ ((2/z) — 1). 

To summarize, we studied in this paper the magnetic phases of the 2D Hubbard model 
at low doping. We found that the planar spiral (ir, Q) phase has lower energy than the Neel 
and (Q, Q) phases in some range of a phenomenological parameter z. However, this state 
has a negative static susceptibility in a region of q space around (tt, 7r), and is therefore 
unstable. We searched for the stable state, and found that it corresponds to a very different 
noncoplanar spin configuration. This, in turn, gives rise to a novel scenario of the trans- 
formation of the equilibrium magnetic state with doping. Namely, upon doping, the initial 
commensurate antiferromagnetic ordering remains unchanged except for a small decrease in 
the amplitude of the order parameter. At the same time, doping gives rise to a transverse 
component of the order parameter which forms a spiral in the plane perpendicular to the 
direction of commensurate order. This transverse component is small to the extent of x, 
and the low-T behavior at finite doping remains nearly the same as in the commensurate 



antiferromagnet POfl. 



Finally, we discuss how these results can be applied to the experimentally studied disor- 
dered regions of cuprate superconductors. The hope here is that the strongest divergence in 
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the susceptibility in the ordered state corresponds to the maximum value of the peak in x 
in the region where the correlation length is finite. In the absence of a long-range order, the 
actual susceptibility is a mixture of x + ~ and X zz ■ In the spiral state, both x H and x zz had 
peaks at the incommensurate momenta ±Q. In the noncoplanar state however, the peak 
in x at the commensurate (tt, 7r) point is much stronger than at the incommensurate mo- 
menta. We thus observe that although incommensurability in the nonplanar state is much 
larger than in the spiral state (inverse pitch of the spiral scales as yfx), experimentally it 
is very difficult to distinguish the noncoplanar state from the Neel state. The noncoplanar 
state therefore should yield the same T— dependence of spin-lattice relaxation rate as the 
(jr, 7r) phase. This is consistent with the low-frequency NMR experiments ||. The situation 
at higher frequences (where neutron-scattering experiments were done 0) is less clear as 
high-frequency bosonic excitations in the spiral phase are positive, and this phase may be 
advantageous over noncoplanar state. The possibility of the frequency crossover between 
the two phases is now under study. Note in passing that for U ~ Qt relevant to experiments, 
we obtain q m 6x for the planar spiral state in a striking agreement with neutron scattering 
experiments on La 2 - x Sr x Cu04 J7| which yield q/n « 2x. 
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FIGURES 

FIG. 1. a) Spin configuration of a noncoplanar state. Arrows with thick ends point out of the 
plane, while those with thick tails — into the plane. This configuration is different from the double 
spiral considered in ||b. b) Two adjacent spins in the equilibrium configuration. The in-plane 
component, S± ~ x 1 / 2 , is small compared to the off-plane component, Su. 
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